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ABSTRACT 

A combined analysis of Cosmic Microwave Background (CMB) and Lyman-a forest 
data can constrain the matter power spectrum from small scales of about l/i -1 Mpc 
all the way to the horizon scale. The long lever arm and complementarity provided by 
such an analysis has previously led to a significant tightening of the constraints on the 
shape and the amplitude of the power spectrum of primordial density fluctuations. 
We present here a combined analysis of the WMAP three year results with Lyman- 
a forest data. The amplitude of the matter power spectrum og and the spectral index 
ris inferred fro m the joint analy sis of high and low resolution Lyman-a forest data 
as analysed bv IViel fc Haehnelt (2006) are consistent with the new WMAP results 
to within la. The join t analysis with the mainly low resolution data as analysed by 
McDonald et al. ( 2005) suggest a value of erg that is ~ 2a higher than that inferred 
from the WMAP three year data alone. The joint analysis of the three year WMAP 
and the Lyman-a forest data also does not favour a running of the spectral index. 
The best fit values for a combined analysis of the three year WMAP data, other CMB 
data, 2dF and the Lyman-cv forest data are (erg, n s ) = (0.78 ± 0.03, 0.96 ± 0.01). 

Key words: Cosmology: observations - cosmology: theory - cosmic microwave back- 
ground, cosmological parameters - quasars: absorption lines 
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1 INTRODUCTION 

Measurements of the matter power spectrum from Lyman- 
a forest data extend to smaller scales and probe a redshift 
range complementary to estimates of the matter power spec- 
trum from Cosmic Microwave Background (CMB), galaxy 
surveys or weak gravitational lensing o bservations (e.g. 
ICroft et all Jl99Sft._lGnedin fc Huil lll99SlhJMcDonald et al] 
feOOCth iHui et al. " 
iMcDonaldl J2003T) 
]2004h . lVielet al.N2004 VHS)). 

The combined analy sis of Lyman-q fores t data with the 
first year of WMAP data JSpergel et alJ200l WMAP1) sug- 
gested that the fluctuation amplitude of the matter power 
spectrum on small scales was rather high (as ~ 0.9) and 
that there was no significant deviation of the spectral in- 
dex of primordial density fluctuations from a Harrison- 
Zeldovich spectrum (ris = 1). There was also no evidence 
for a (large) running of the spectral index, a non-zero neu- 
trino mass or a devia tion from a cold dark m atter spec- 
trum a t small scales J Viel et al] J20041 VHS). IViel et al] 
l2004). iMcDonald et al] J2005I. M0 5 ). IViel et al] fc OOE tT 

[ Seliak et al.l ll2005lhlLidz et alJlbooilhlBeltran et all 1200*3) . 
Abazaiianl l200d) h 

IViel et al] J2004D found as = 0.94±0.08, n = 0.99±0.03 



(la) and no evidence for a (large) running of the spectral in a 
combined analysis of a large sample of high resolutio n spec- 
tra quasar (QSO) abso rption spectra at z ~ 2.5 (Ki m et al] 
12004: ICroft et al]l2002l) and the WMAP1 data. Similar re- 
sults, with somewhat smaller errors (erg = 0.90 ± 0.03, 
ris = 0.98 ± 0.02) have been subsequently obtained by the 
SDSS collaboration in a combined analysis of the WMAP1 
and other CMB data , SDSS galaxy sur vey data and SDSS 
Lyman-a forest data JSeliak e t alj[20 05T). The Lym an-q for- 
est data analysed by M05 and ISeliak et al] J2005I) consists 
mainly of low-resolution low S/N SDSS spectra with a wide 
redshift coverage (2 < z < 4) to which they ad ded a small 
samp le of eight high resolution spectra (McDonald ct al. 
2000). The flux power spectrum was modelled using dark 
matter simulations, which take into account hydrodynami- 
cal effects in an approximate way and were calibrated with 
a few hydrodynamical simulations. IViel fc Haehnelt! (J2006) 
found a 8 = 0.91±0.07, n s = 0.95±0.04 for the SDSS Lyman- 
q forest data alone using a suite of state-of-the-art full hy- 
drodynami cal simulations. Further stud i es of Lyman - q for - 
est data bvlDesiacques fc Nusseil ]2005l) . Ijena et al] (12005]) 
and lZaroubi et al] (l2005fl also came to similar conclusions. 

The WMAP3 data alone argues now for significant devi- 
ation from a Harrison-Zeldovich spectrum, ris = 0.95 ± 0.02, 
and a smaller value for the fluctuation amplitude on small 
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Figure 1. 1 and 2<r likelihoods for eg and n^ marginalized over all other parameters. Left panel: Constraints are for WMAP1 only 
(green), the LUQAS+CROFT data sets as analysed by VHS (cyan) and the SDSS Lyman-a forest data of M05 (blue). The thick dashed 
white empty contours refer to WMAP1 + VHS, while the solid blue contours are for WMAP1 + SDSS. Right panel: As in the left panel 
but for the WMAP3 data set. 



scales (T 8 = 0.74 ± 0.06 ijSnergel et al.ll200rt . The WMAP 
team chose not to update their combined analysis of CMB 
and Lyman-a forest in their WMAP3 data release. In this 
Letter we will present such a joint analysis. 



2 THE DATA SETS 

2.1 WMAP 

The WMAP 1 satellite has mapped the entire sky in five 
frequency bands between 23 and 94 GHz with polarization 
sensitive radiometers. The temperature power spectrum has 
been measured over a lar ge range of scales (I < 1000) to an 
unprecedented accuracy ( Hi nshaw et alj|2 006l) . We will use 
the temperature and polarization dPage et alJl2006l) power 
spectra and maps as used by the WM AP likelihood codes 
i Verde et alj 120031 : ISpergel et alJl2006| ) as im plemented in 
the code COSMOMC 2 JLewis fc Bridlell2002l) . 



2.2 Lyman-a forest data sets 

We will investigate two different Lyman-a forest data sets. 
The sample o f high resol u tion QSO absorption spectra used 
by VHS and IViel et alJ J2004I) . c onsisting of the LUQAS 
sample (27 high resolution QSOs) llKim et al.ll2004l) and the 
(reanalysed) sample of C02 (30 high resolution and 23 low 
resolution spectra), and the SDSS Lyman-a forest data as 



presented bv iMcDonald et al.1 (120061) . The SDSS Lyman- 
a forest data set consists of 3035 QSO spectra with low 
resolution (R ~ 2000) and low S/N (~ 10 per pixel) span- 
ning a wide range of redshifts (z — 2.2 — 4.2), while the 
LUQAS and the C02 samples contain mainly high resolu- 
tion (R ~ 45000), high signal-to-noise (> 50 per pixel) QSO 
spectra with median redshifts of z — 2.125 and z — 2.72, re- 
spectively. Modelling the flux power spectrum of the Lyman- 
a forest accurately for given cosmological parameters is not 
as straightforward as modelling the CMB power spectra 
and accurate numerical simulations are required. M05 mod- 
elled the flux power spectrum using a large number of Hy- 
dro Particle Mes h (HPM) simulations dGnedin fc Huilll998l : 
IViel et alJ 2006), calibrated with a few full hydrodynam- 
ical simulations. VHS improve d instead the effective bias 
met hod developed by C02 (s ee lGnedin fc Hamilton! J2002ft 
and IZaldarriaea et a l. (2003) for a critical assessment of 
the errors involved), by using a grid of full hydrodynam- 
ic al simulations run with the Tree -SPH code GADGET- 
2 ijSpringel et al.l200ltlSp ringcl 2005) to infer the linear mat- 
ter power spectrum. IViel fc Haehneltl (120061) used a Taylor 
expansion of the flux power spectrum around best fitting 
values based on full hydrodynamical simulations to model 
the dependence of the flux power on cosmological and as- 
trophysical parameters in their independent analysis of the 
SDSS Lyman-a forest data. 



1 http://lambda.gfsc.nasa.gov 
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Figure 2. 1- and 2a likelihoods for ag and Us marginalized over all other parameters. Left panel: The blue contours show the constraints 
for SDSS only as analysed by M05 wit h an HST prior 72 ± 8 km/s/Mpc (blue). The empty solid contours are for the SDSS data set 
as analysed by Vicl & Haclmclt (2006, SDSS-d). The WMAP3 results are shown in green. Right panel: The combined analysis with 
WMAP3 using the same colour/line coding as in the left panel. 



Table 1. Summary of the constraints on ag and ns for the V HS IViel et al.ll2004) and SDSS JMcDonald et alJl200fi> samples. SDSS-d 
refers to the analysis of SDSS data bv lVie l & Haclmclt 12006}). The quoted errors are the 68% confidence limits. 



WMAPl+VHS WMAP1 + SDSS 



WMAP3 + VHS WMAP3 + SDSS WMAP3 + SDSS-d 



er 8 0.90 ±0.10 
ns 0.99 ± 0.04 



0.94 ± 0.08 
0.99 ± 0.03 



0.93 ±0.04 
0.99 ±0.03 



0.76 ±0.05 
0.96 ±0.02 



0.78 ±0.05 
0.96 ± 0.02 



0.86 ± 0.03 
0.96 ± 0.02 



0.80 ±0.04 
0.96 ±0.01 



3 RESULTS 

3.1 Incorporating the Lyman-a data into 
COSMOMC 

The linear dark matter power spectra inferred from the 
two Lyman-a forest data sets have been incorporated into 
the new public available version of COSMOMC. The VHS 
power spectrum consists of estimates of the linear dark mat- 
ter power spectrum at nine values of the wavenumber k at 
z — 2.125 and nine values at z — 2.72, in the range 0.003 < k 
(s/km)< 0.03. The estimate of the uncertainty of the overall 
amplitude of the matter power spectrum is 29%. This esti- 
mate takes into account possible systematic and statistical 
errors (see the relevant tables of VHS for a detailed discus- 
sion). M05 provided a measurement of slope and amplitude 
of the matter power spectrum at z — 3 at a wavenumber 
k = 0.009 s/km with an estimate of the la error of the am- 
plitude of ~ 14%. M05 have also made available a table that 
gives the minimum \ 2 for a given cosmological model as a 
function of the amplitude and slope after marginalization 
over a wide range of cosmological and "nuisance" param- 
eters. The nuisance parameters characterize a range of as- 
trophysical and noise-related systematic uncertainties. We 



have furthermore implemented the model ling of the SDSS 
flux p ower spectrum with the method of IViel fc Haehneltl 
(2006) based on a Taylor expansion of the flux power spec- 
trum around a best fitting model. 



3.2 Constraints on ag and 74 

To make contact with previous analyses we show the 
marginalized la and 2a likelihoods in the <rg-n s projection 
in Figure 1. The coloured contours in the left panel show 
the constraints for the VHS sample (light cyan), the SDSS 
sample (dark blue) with the likelihood estimates provided 
by M05 and the constraints for the WMAP1 data (light 
green). In all cases we assume the Universe to be flat, no 
contribution from tensor perturbations, a pure cosmological 
constant (w = — 1) and neutrinos with negligible mass. For 
the analysis of the WMAP1 data we assumed a prior on the 
Thomson optical depth < r < 0.3. 

The dark solid (SDSS) and light dashed (VHS) con- 
tours show the joint constraints fo r the Lyman-a forest an d 
WMAP1 data. As pointed out bv lViel fc HaehnebJ feOOfJ) . 
there is remarkable agreement between the two joint anal- 
yses with the high resolution absorption spectra as anal- 
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ysed by VHS and I Viel et al.l <I2004I) and the larger sample 
of low resolution SPSS spectra as analysed by M05 and 
ISeliak et al.l 1)20051) . The Lyman-a forest data break some 
of the degeneracies of the WMAP1 data which are responsi- 
ble for the elongated shape of the error contours. The joint 
analysis tightens the constraints in the erg — iin plane by a 
factor ~ 2 (VHS) and ~ 4 (M05), compared to the con- 
straints from the WMAP1 data alone but offers less help 
in improving the constraints on ris. As discussed above and 
summarized in Table 1, in the combined analysis with the 
WMAP1 data the best fitting value of the spectral index is 
not significantly different from n s = f and as ~ 0.9. The sig- 
nificantly smaller error bars of the SDSS data set are due to 
the much larger sample: the wider range of redshifts covered 
is responsible for breaking some of the degeneracies intrinsic 
to the Lyman-a forest data (M05). 

In the right panel of Figure Q we show how the situ- 
ation has changed w ith the WMAP3 data . As discussed in 
ISpereel et alJ {2006) and lPaee et alj (l2006f) . the contours for 
the WMAP3 data have shrunk by a factor ~ 3 compared to 
those of WMAP1 and lie at the bottom part of the region 
allowed by the latter. This is mainly due to the improved 
measurem ent of the optical depth from the large scale po- 
larization dPage et al.ll2006l:lLewisll2006t) . 

In the case of the high resolution VHS sample the 
errors are too large to significantly tighten the WMAP3 
constraints. The joint analysis of WMAP3 with the SDSS 
Lyman-a forest data places now even tighter constraints at 
the bottom end of the range preferred by the SDSS Lyman- 
a forest data alone, about 2a above the best fit valu e from 
the W MAP3 only data: <j 8 = 0.86 ± 0.03 (see also iLewis! 
(2006)). The discrepancy is larger than may be naively in- 
ferred from the overlap of the Lyman-a only and WMAP3 
only analysis because the data sets prefer different values 
for some of the other parameters in particular f2 m . The best 
fit value for fi m (~ 0.3) is also ~ 2a higher for the com- 
bined analysis than for the WMAP 3 data alone. The SPSS 
Lyman-a forest data as analysed bv lMcDonald et al.l (l2005h 
and the new WMAP results appear to be marginally con- 
sistent. Possible explanations for the (small) discrepancy, if 
there exists one, may be somewhat too optimisti c errors for 
one o r both of the data sets. As discussed in IPage et all 
J2006TI the polarization measurements of the CMB are very 
difficult mainly due to foreground polarization. Moreover, 
modelling the Lyman-a forest data also has a range of not 
yet fully understood systematic unc ertainties. 

As discussed exte nsively by JMcDonald et al.l 120051) 
IViel fc Haehnel t (2006) the mayor systematic uncertainties 
are the still not very accurately known thermal state of 
the IGM, the determination of the effective optical depth, 
the m odelling of the e ffect of strong absorption line sys- 
tems (IViel et alJ i2004TO and UV fluctuations and the re- 
maining deficiencies in our ability to accurately predict the 
flux power spectrum for a large parameter space. To inves- 
tigate the last issue issue further we compare in the left 
panel of Fig ure 2 the analysis of th e SDSS data by M05 
with that by IViel fc Haehneltl J2006T) which is based on a 
Taylor expansion of the flux power spectrum around a best 
fit ting model (labell e d as S DSS-d in Figure 2). The analysis 
of IViel fc Haehneltl J2006t) uses more accurate full hydro- 
simulations instead of the approximate simulations of M05 
at the expense of a much less complete sampling of param- 



eter space, especially far from the best-fitting values, where 
the errors are possibly underestimated. Note that here, 
both for the forest data and the joint anal ysis, we assumed 
an H ST prior for the Hubble constant IFreedman et alJ 
200f) that significantly shrinks the error contours for the 
Lyman-a for est data compared to that of Figure f. As 
discussed by IViel fc Haehneltl (I2006f) there is remarkable 
agree ment between the two a nalyses (note that the anal- 
ysis of IViel fc Haehneltl (2006) does not use the last three 
redshift bins nor high resolution QSOs compared to that 
of M05). In the right panel of Figure 2 we show the con- 
straints for the joint analysis with the W MAP3 data. The 
SDSS Lyman-a forest data as analysed bv lViel fc Haehneltl 
(2006) combined with WMAP3 give a smaller value best 
fit value of as = 0.80 ± 0.04 which is in agreement with 
that from the WMAP3 data alone to within la. The joint 
analysis of Lyman-a forest data and the new WMAP data 
including a possible running of the spectral index gives 
n run = -0.002 ±0.01 5 at fc = 0.002 Mpc" 1 also in agreement 
with the estimate bv lSpergel et al.l (2006). 

We have also performed an extended combined 
analysis that includes t he furt her CMB experiment s 
ACB AR <Kuo et al.l \2004) CBI (iReadhead et alJ l2004Tl . 
VSA {Dickinson^^y J2OO4J) , the 2dF galaxy power spec- 
trum dPercival et all 120011) and the VHS and SDSS-d 
Lyman-a forest data. In this case we get (as, rh) = (0.78 ± 
0.03, 0.96 ± 0.01 ). Further results are listed in Tabled 



4 CONCLUSIONS 

We have performed a combined analysis of the WMAP three 
year results with high and low resolution Lyman-a forest 
data in order to constrain the shape of the power spectrum 
of primordial density fluctuations and the amplitude of the 
matter power spectrum at intermediate scales as- The main 
results are as follows. 

• The high resolution VHS Lyman-a forest data is con- 
sistent to within la with the three year WMAP results but 
offers little additional constraining power due to the large er- 
ror bars. The larger sample of mainl y low resolution L yman- 
a forest data (SDSS) as analysed bv lViel fc Haehn elt (2006) 
is also consistent to within la with the new WMAP results. 
However the joint analysis of the SDSS data gives about 2a 
higher as and f2 m values than those inferred from the new 
WMAP results alone. 

• The best fit values for a combined analysis are (la): 
(a 8 , Th) = (0.78 ± 0.05, 0.96 ± 0.02) and (as, n s ) = (0.86 ± 
0.03, 0.96 ± 0.02) for WMAP combined with high resolution 
Lyman-a forest data and WMAP combined with low reso- 
lution Lyman-a forest data as analysed b y M05. The anal- 
ysisof the SDSS data set as analysed bv lViel fc Haehneltl 
(2006) based on full hydrodynamical simulations gives as = 
0.80 ± 0.04. 

• The joint analysis of Lyman-a forest data and the new 
WMAP data does not favour a running of the spectral index. 
The best fitting value is n run = -0.002 ± 0.015 at fc = 0.002 
Mpc" 1 . 

• Adding other CMB data sets, the 2dF galaxy survey 
and both Lyman-a data sets the constraints on the matter 
power spectrum become (as, ns) — (0.78±0.03, 0.96±0.0f). 
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Table 2. The marginalized constraints on cosmological parameters from WMAP3 and other data sets. VHS refers to the LUQAS+CROFT 
sample as analysed b^J^iel^^al. (2004); SDSS refer s to the mea sureme nt bvlM cDonald ct al. (2005); EXT refers to small er scales CMB 
data sets: ACBAR JKuo et all2004l) . CBI faeadhead et all2004l). VSA JDickinson et all2004l) and the 2dF galaxy survey iPercival et alj 
120011) . SDSS-d refers to the SDSS analysis bv ) Viel fc Haehneltj J200rj) . 



WMAP3+VHS 



WMAP3 + SDSS 



WMAP3 + SDSS-d WMAP3 + EXT + SDSS-d WMAP3 + SDSS-d(ruiin.) 



n c h 2 


0.106 ±0.007 


0.109 ±0.008 


0.120 ±0.006 


0.110 ±0.006 


0.109 ± 0.006 


0.112 ±0.007 


io 2 n h h 2 


2.222 ± 0.069 


2.237 ±0.072 


2.277 ±0.065 


2.226 ±0.071 


2.224 ±0.066 


2.221 ± 0.094 


1 £ m 


0.242 ± 0.032 


0.257 ±0.037 


0.304 ±0.031 


0.258 ±0.029 


0.253 ± 0.028 


0.269 ±0.036 


h 


0.729 ±0.029 


0.719 ±0.031 


0.688 ± 0.025 


0.719 ±0.026 


0.723 ± 0.025 


0.711 ±0.034 


T 


0.089 ± 0.030 


0.092 ± 0.029 


0.101 ±0.028 


0.098 ± 0.032 


0.104 ±0.036 


0.104 ±0.030 


CT8 


0.761 ±0.046 


0.784 ± 0.048 


0.857 ±0.028 


0.801 ±0.039 


0.785 ± 0.035 


0.800 ± 0.037 


«s 


0.956 ±0.016 


0.956 ±0.017 


0.964 ± 0.016 


0.960 ± 0.013 


0.957 ±0.014 


0.963 ± 0.020 


^run 




- 


- 


- 


- 


-0.002 ± 0.015 



The Lyman-Q forest data appears to be in reasonable 
agreement with the CMB and other data sets which probe 
the matter power spectrum at larger scales. The Lyman- 
q forest data will thus continue to unfold its special power 
to measure parameters that affect the overall shape and/or 
the small scale part of the matter power spectrum. For the 
near future further progress is, however, likely to be driven 
by a better understanding of the systematic uncertainties, 
rather than the compilation of larger data sets. 
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